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Combined Interference Cancellation and Maximum 
Likelihood Decoding of Space-Time Block Codes 

5 

Reference to Related Applications 

This application claims the benefit of U.S. Provisional Application 
No. 60/061,145, filed October 6, 1997. This application is also related to U.S. 
Application 09/074224, filed May 7, 1998, titled 'Transmitter Diversity 
10 Technique for Wireless Communications". 

Background of the Invention 

This invention relates to wireless communication and, more 
particularly, to techniques for effective wireless communication in the 

15 presence of fading, co-channel interference, and other degradations. 

Rapid growth in mobile computing and other wireless data services is 
inspiring many proposals for high speed data services in the range of 64-144 
kbps for micro cellular wide area and high mobility applications and up to 2 v 
Mbps for indoor applications. Research challenges include the development 

20 of efficient coding and modulation, and signal processing techniques to 
improve the quality and spectral efficiency of wireless communications and 
better techniques for sharing the limited spectrum among different high 
capacity users. 

The physical limitation of the wireless channel presents a 
25 fundamental technical challenge for reliable communications. The channel 
is susceptible to time-varying noise, interference, and multipaths. Power and 
size limitations of the communications and computing device in a mobile 
handset constitute another major design consideration. Most personal 
communications and wireless services portables are meant to be carried in a 
30 briefcase and/or pocket and must, therefore, be small and lightweight. This 
translates to a low power requirement since small batteries must be used. 
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However, many of the signal processing techniques which may be used for 
reliable communications and efficient spectral utilization demand significant 
processing power, precluding the use of low power devices. Continuing 
advances in VLSI and integrated circuit technology for low power 
5 applications will provide a partial solution to this problem. Still, placing 
most of the signal processing burden on fixed locations (base stations) with 
relatively larger power resources than the mobile units will, likely, continue 
to be the trend in wireless systems design. 

Perhaps the single most important parameter in providing reliable 
TO communications over wireless channels is diversity. Diversity techniques 
which may be used include time, frequency, and space diversity 

• Time diversity: channel coding in combination with limited interleaving is 
used to provide time diversity. However, while channel coding is extremely 
effective in fast fading environments (high mobility), it offers very little 

15 protection under slow fading (low mobility) unless significant interleaving 
delays can be tolerated. 

• Frequency diversity: the fact that signals transmitted over different 
frequencies induce different multipath structures and independent fading is. 
However, when the multipath delay spread is small compared to the symbol 

20 period, frequency diversity is not helpful. 

• Space diversity: the receiver/transmitter uses multiple antennas that are 
separated for reception/transmission and/or differently polarized antennas to 
create independent fading channels. Currently, multiple antennas at base- 
stations are used for receive diversity at the base. However, it is difficult to 

25 have more than two antennas at the mobile unit due to the size and cost of 

multiple chains of RF conversion. 

Previous work on transmit diversity can be classified into three broad 

categories: schemes using feedback, schemes with feedforward or training 

information but no feedback, and blind schemes. The third category (blind 
30 schemes) relies primarily on multiple transmit antennas combined with 

channel coding to provide diversity. An example of this approach is 
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disclosed in the aforementioned copending application 09/074224, 1998, 
titled "Transmitter Diversity Technique for Wireless Communications," filed 
May 7. 

5 Summary of the Invention 

Improved performance is attained in an illustrative arrangement 
where K synchronized terminal units that transmit on N antennas to a base 
station having M > K antennas, by combining interference cancellation (IC) 
and maximum likelihood (ML) decoding. More specifically, space-time 
10 block coding is employed in transmitters that employ N transmit antennas 
each, and the signals are received in a receiver that employs M receiving 
antennas. In accordance with the processing disclosed herein, by exploiting 
the structure of the space-time block code, K-l interfering transmitting units 
are cancelled at the receiver, regardless of the number of transmitting 
15 antennas, N, when decoding the signals transmitted by a given mobile unit. 
In another embodiment of the principles of this invention, signals of a first 
terminal unit are decoded first, and the resulting decoded signals are 
employed to cancel their contribution to the signals received at the base t 
station antennas while decoding the signals of the remaining AM terminal 
20 units. The process is repeated among the remaining K-l terminal units. That 
is, among the remaining AM, signals of a first terminal unit is decoded first 
and the resulting decoded signals are employed to cancel their contribution to 
the signals received at the base station antennas while decoding the signals of 
the remaining K-2 terminal units, and so on. This procedure is repeated M 
25 times, each time starting with decoding signals of a particular terminal unit. 
This successive procedure will yield additional performance improvement. 

Both zero-forcing (ZF) and minimum mean-squared error (MMSE) 
interference cancellation (IC) and maximum likelihood (ML) techniques are 
disclosed. 



30 



Brief Description of the Drawing 
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FIG. 1 depicts an arrangement that, illustratively, includes a receiving 
base station (20) and two transmitting terminal units (10 and 30), 

Detailed Description 

5 FIG. 1 illustrates two transmitting units and one receiving unit that 

comport with the principles disclosed herein. However, this is merely 
illustrative, and the disclosed method is useful for more than two terminals 
{K>2\ 

10 Single Transmitting Unit 

Transmitting unit 10 may correspond to the transmitting circuitry in a 
terminal unit, while receiving unit 20 may correspond to the receiving 
circuitry in a base station. Terminal unit 30 is shown identical to terminal 
unit 10. It should be understood, of course, that each terminal unit has a 

15 receiving circuit, and the base station has a transmitting circuit. The terminal 
units are shown to have two antennas each. Receiving unit 20 is also shown 
to have two receiving antennas. Here, too, it should be kept in mind that, 
generally, any number, M > 2, of receiving antennas can be had. Particular 
advantage is realized when M> K. Since the mathematical treatment below 

20 is couched in general matrix notations, the expressions are valid for any 
number K and/or M 

Considering terminal unit 10, the information source provides input 
symbols to element 13 which develops a block code. The symbols are 
divided into groups of two symbols each, and at a given symbol period, the 

25 two symbols in each group {c l9 c 2 } are transmitted simultaneously from the 
two antennas. The signal transmitted from antenna 1 1 is c, and the signal 
transmitted from antenna 12 is c 2 . In the next symbol period, the signal 
- c, * is transmitted from antenna 1 1 and the signal c, * is transmitted from 
antenna 12. The symbols are modulated prior to transmission with 

30 constellation mappers 14 and 15, followed by pulse shapers .16 and 17, 
respectively, in a conventional manner. 
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In receiver 20, signals are received by antennas 21 and 22 and are 
applied to detector 25. 

In the mathematical development of the algorithms disclosed herein, 
it is assumed that the channel from each of the two transmit antennas remains 
5 fixed over two consecutive symbol periods. That is, 

h s {nT) = />,((» + 1)D, ■/ = O) 
To ascertain the channel characteristics, the transmitter carries out a 
calibration session, during which pilot signals or tones are transmitted. The 
signals received during the calibration session are applied to channel 
10 estimator circuits 23 and 24, which are well known circuits, and the channel 
characteristics are thus obtained. - 

When only transmitter 10 is considered, the received signals at 
antenna 2 1 can be expressed as 

r M = h n c } +h l2 c 2 +r), (2) 
15 r l2 =-h u c 2 + h l2 c x +r} 2 (3) 

where /'// and r\i are the received signals over two consecutive symbol 
periods, h u denotes the fading channel between transmit antenna 11 and 
receive antenna 21, h xz denotes channel between transmit antenna 12 and 
receive antenna 21. and tj x and tj 2 are noise terms, which are assumed to be 
20 complex Gaussian random variables with zero mean and power spectral 
density No/2 per dimension. Defining the vectors r = [r xl r l2 *] T , c = [c,c 2 ] r , 
and 77 ~ [tj x ti 2 *] t , equations (2) and (3) can be rewritten in a matrix form as 

r = H-c + 7 (4) 
where the channel matrix H is defined as 



25 H = 



h * -h * 



The vector T| is a complex Gaussian random vector with zero mean 
and covariance N 0 > I . Defining C as the set of all possible symbol-pairs 
c = {c, 3 c 2 } , and assuming that all symbol pairs are equi-probable, it can be 



(5) 
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easily shown that the optimum maximum likelihood (ML) decoder selects 
from C the symbol-pair c that minimizes the expression ||r — H-cJ| 3 . This 
can be written as 

c = arg mm||r-H*c||\ ( 6 ) 

It was shown by S. Alamouti in "Space Block Coding: A simple 
Transmitter Diversity Scheme for Wireless Communications/* submitted to 
IEEE JSAC, September 1 997 that the diversity order of the above space-time 
block code is equivalent to thai of a two branch maximal ratio receive 
combining (MRRC). Because of the orthogonality of the matrix H, 
Alamouti also showed that this decoding rule decomposed into two separate 
decoding rules for C\ and ci. The uncertainty, A r , of the decoded symbols i 
is defined as 



A. = r-H-c 



, and H 2 = 



h ' * — h * 



(7) 



15 It should be noted that the above analysis addresses only receive 

antenna 21. When receiver 20 uses both antennas, i.e., antennas 21 and 22, 
two received signal vectors rj and r? can be defined for antenna 21 and 22, 
respectively, as 

r ( =H, -c + 7, (8 ) 
r 2 = H 2 *c + 7 2 ^ 
where H| and H : are the channel matrices to receive antennas 21 and 22, 
respectively, and 77, and rj 2 are the corresponding noise vectors. That is, 



(9a) 



where h 2l denotes the channel between transmit antenna 12 and receive 
25 antenna 22, and h 2Z denotes the channel between transmit antenna 1 1 and 
receive antenna 22. In this case, the ML decoding rule is 

c = arg mintjjr, -'H, *c|| 2 -f jjr 2 - H 2 • c|| 2 ) , ( 1 0) 
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A^ih-H^cH^lir^H^cI 2 . 



(11) 



10 



15 



20 



As before, both the matrices H, and H 2 are orthogonal matrices and hence 
the above decoding rule also decomposes to two separate decoding rules for 
Cj and c 2 . Note that the rate of transmission (of information symbols) in the 
space-time block coding scheme is 1 

Interference Cancellation and ML Decoding: BASIC CASE 
FIG. 1, however, shows two terminal units, and the issue that needs 
to be addressed is the detection performance at the base station receiver 
when the two terminal units transmit over the same time and frequency 
channel. 

In the notation below, g u denotes the fading channel between 
transmit antenna 31 and receive antenna 21, g l2 denotes the channel 
between antenna 31 and antenna 21, g 2l denotes the channel between 
antenna 32 and antenna 22. and g 22 denotes the channel between antenna 32 
and antenna 22. Also, {c r c 2 } and {s,,s 2 ] denote the two symbols 
transmitted from terminal units 10 and 30, respectively. 

At receiver 20, the received signals over two consecutive symbol 
periods at receive antenna 21, r n and r 12 , are 



Defining r, =[r„r I2 *] r , c = [c, c 2 ] r , s = [s, s 2 ] T \ and n, =[t| n r\ l2 ] T 
equations (12) and (13) can be rewritten in matrix form as 



where the channel matrices H, and G, between the transmitter units 10 and 
30 and receive antenna 2 1 are given by 



r X2 = ~V C 2 +/l 12 C » ~gl\ S l +Sl2 5 I + ^12 



(13) 



(12) 



r. - H, -c + G. *s + n 



(14) 




(15) 
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The vector n, is a complex Gaussian random vector with zero mean and 
covariance N 0 * I . Similarly, the received signals over two consecutive 
symbol periods at receive antenna 22, / 2 , and /%, are 

r :i = h 2x c x + h zz c 2 +#2,5! + g 22 s 2 + ti 2j 
r 22 = -h 2X c[ + h 22 c\ - g 2] s* 2 + g 22 s," + *n 22 

In a similar fashion, defining r 2 = [r 2l r 22 *] T and n : =[J7 2 i , 722 *3 r equations 
(16) and (17) can be rewritten as 

r 2 = H 2 c + G 2 s + n 2 
where the channel matrices H : and G : berween transmitter units 10 and 30 
and antenna 22 are given by 



h,-[.V 



, and G 2 = 


&2 1 &22 
* • 




.§22 *~<§21 _ 



Equations (14) and (18) can be combined to yield the matrix form 

n, 



[rj [h, G,][s_ 



n, 



(16) 
(17) 



(18) 



(19) 



(20) 



15 Zero-Forcing IC and ML Decoding Scheme: In this case, a matrix W 

is chosen such that 



W 





H o" 




c 








0 G 




s 


+ 


- 5 2- 



(21) 



20 



We can find an appropriate matrix W by realizing that 



th, G.rrAr 1 oir i -c.g;'] 

|H, GJ [ 0 Aj'JL-HjHr 1 I J 



where 

A, =H, -G,G;'H : and A, = G : - H_,H7'G, 
Hence, if we select Was 



i -g,g; 
-h,h;' i 



(22) 



(23) 



(24) 
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we will have 



W 



;]-R]-[f i. 



n, 



(25) 



10 



15 



20 



where 

G = G 2 — H 2 Hj 'Gj 

iij = Dj — GjGj'n., 

n., = n : — H 2 H, 'n, 
From equation (25) it can be easily observed that the modified 
received signal vector r ] contains signals only from transmitter 10 ( i.e. 
signals from transmitter 30 have been canceled or removed) and, 
correspondingly, the modified received signal vector ? : contains signals only 
from transmitter 30 ( i.e. signals from transmitter 10 have been canceled or 
removed). A number of other attributes can be shown from the above to be 
true. 

1) the modified noise vector n : is a zero mean complex Gaussian random 
vector with covariance 



i + - 



D m 



I 



^2 J 



where D gi = |g n | 2 + |g 12 | 2 and D g2 = \g 2l \ 2 +|^ 22 | 2 . Hence, the modified 
noise vector n, is also white. 

2) the modified noise vector n 2 is also a zero mean Gaussian random vector 
with covriance 



D, 



I 



A I J 



where D M = |/7,,| 2 + j/'| 2 j" and A: = \ l7 u\~ l^'aa i~ » an< * hence it is also white. 
3) The matrices H and G the form 



H = 



h' -JL' 



, and G = 



g; g 2 
£2 ~g> J 



(26) 



(27) 



(28) 



(29) 
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4) Conditioned on G, and G : , the random variables h x and h 2 are both zero 
mean complex Gaussian random variables with variance a* = 1 + ^ S \/D g2 . 

5) Conditioned on H, and H,, the random variables g x and g 2 are both zero 
mean complex Gaussian random variables with variance = 1 + D h2 /D kx . 

5 6) The modified channel matrices H and G have a structure similar to that in 
equation (5), i.e. the modified channel matrices H and G are orthogonal 
matrices. 

Considering the modified received signal vector r t , which contains 
signals only from transmitter 10, i.e., 
10 r; = H - c + n,, {30) 

it is noted that this expression resembles the expression in equation (4). 
Hence, in accordance with the principles disclosed therein, the optimum ML 
decoder for the symbols transmitted by terminal unit 1 0 evaluates an 
equation that is similar to the expression of equation (6), and is given by 

15 c = arg min||if - H • cf . (3 ! ) 



The corresponding decoder uncertainty is given by 



Moreover, since the channel Matrix H is orthogonal, the ML decoder will 
also decompose into two separate rules for c\ and c 2 . 
20 In a similar fashion, considering the modified received signal vector 

r 3 , which contains signals only from transmitter 10, i.e., 
K = H c + H 2 , 

it is noted that this expression resembles the expression in equation (4). 
Hence, the optimum ML decoder for the symbols transmitted by terminal 
25 unit 30 evaluates an equation that is similar to the expression of equation (6), 
and is gi ven by 

s = argmin te - G • ill . 
The corresponding decoder uncertainty is given by 



(32) 



(33) 



(34) 
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A s =||r 2 -G-s|| 2 .' (35) 

Moreover, since the channel Matrix G is also orthogonal, the ML decoder 
will also decompose into two separate rules for S\ and s 2 . 

The above-disclosed technique can be easily implemented within a 
detector 25 that comprises a stored program general purpose processor. 
Specifically, a subroutine ( c, A )=ZF.DECODE(r l5 r 2 , H,, H 2 , Gj, G 2 ) can be 
installed which returns the values c,A in response to submitted inputs 
ri,r 2 ,HiJH2,Gi, and G 2s as shown below: 

(c,A,) = ZF.pECODE^pr^HpH^GpG,) _ ; 

t . .. .. - ■ . . - ... 

7 = r 1 -G l q; , r a 

fi-H, -GjG^Hj 
c = are min j| r - H • cf 2 



A r = |r-H.c|| 

} 

With such a subroutine, both s and c can be estimated, as follows: 

(c,A)=ZRDEeODE(r 1 ,r 2 ,H,,H 2 ,G l ,G 2 ) (36) 
(s,A)=ZF.DECODE(r 2 ,r w G 2 ,G I ,H 2 ,H I ). (37) 

15 It may be noted that knowledge of the noise power is not required. 

Simulation results reveal that the performance of the FIG. 1 system which 
employs the principles disclosed herein and the ZF.DECODE subroutine is 
equivalent to that when onty one terminal unit exists and the base station 
uses a single receive antenna. That is, the reception is equivalent to the 

20 performance of two branch MRRC diversity receiver. Advantageously, 
however, the disclosed technique is able to support two co-channel terminal 
units. 

The discussion above centered on describing the technique for 
canceling out the signal of transmitter 10 when detecting the signal of 
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transmitter 30, and for canceling out the signal of transmitter 30 when 
detecting the signal of transmitter 10. Effectively, detector 25 of receiver 20 
can either comprise two processors, with one making the subroutine call of 
equation (31) and the other making the subroutine call of equation (32). 
Alternatively, the signals can be stored within detector 25 and the subroutine 
calls of equations 3 1 and 32 can be made seriatim. 

Minimum Mean-Squared Error IC and ML Decoding Scheme : The 
above-disclosed approach for canceling the contribution of an interfering 
terminal unit is known as the zero-forcing (ZF) as a minimum mean-squared 
error technique (MMSE). 

Recalling equation (20), the vector r can also be written as 
r = H- c + n 



H = [ H ' G '1 



To simplify notations, the vector r is also defined as r = [r, r, r 3 r 4 ] T . 

When seeking to detect and decode signals {c,,c 3 } by minimizing a 
mean-squared error criterion, the goal is find a linear combination of the 
received signals such that the mean-squared error in detecting the signals 
{c,,c 2 } is minimized. In general terms, this can be expressed by an error 
cost function that is to be minimized, such as the function 

J ( a >& = ||i>.^ (P,c, + P 2 c 2 )f = \a ■ r - p . cf. 

II '"I || 

One may note that a minimum is certainly reached when both a and fi are 
equal to 0, but that, of course, is not desired. Therefore, either p, or p 2 is set 
to 1. When p 2 is set to 1, we get the following minimization criterion from 
equation (40) 



wherea l =[a M ,a 12 ,a 13 ,a 14 ,-p 1 ] = [«, _pj and i; =[r 7 ctf . From this it 



(38) 



(39) 



(40) 



(41) 
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can be seen that 



10 



15 



20 





H Q r Tc" 




n 




,o i JU. 


+ 


0_ 



where 0 = [0 0 0- 0]. 

What is needed is to select a x so that the expected value of the 
expression in equation (41) is minimized. That is, select a x to minimize 

E{ J, (5 , )} = E{(a x % - c 2 Xafi - c 2 ) *}. 
Taking the partial derivative with respect to 5, and setting it to zero, what 
results is 



"M h," 










"P.. 




_0_ 



where M = HH* ^ ' s * e signal to noise ratio, I is the 4 by 4 identity 

matrix, h, is the first column of H, and h, is the second column of H. It 
follows that 



It can be shown that 



which yields 



a] =(M-h,h;)" , h 2 and fi\ = h;(M - h,h;)"'h,. 
M-'h.hjM"' 



(M-h.h;) - ' =M -1 — 



l-hjM-'h, ' 



b:M"'h, 



25 



l-h;M"h, 

From the structure of the matrix H we can easily verify that hi and h2 are 
orthogonal. Using this fact and the structure of the matrix M, it can be 
shown that 

P, =0 

a\ = M"'h 2 

Hence, the MMSE IC solution given in equations (45) and (46) will 
minimize the mean-squared error in c, without any regard to c, . Considering 
the alternative cost function when Pi is set to 1, a similar analysis leads to the 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 
(49) 
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conclusion that 

P 2 =0 (47) 
a\ = M-'h, (48) 
The value of V and the values of h.. and g tJ , and consequently, the values of 
H and M are obtained from a training sequence in a conventional manner by 
elements 23 and 24. Since, as indicated earlier, this is quite conventional and 
does not form a part of this invention, for sake of conciseness additional 
details are not presented. In this case, the MMSE IC solution given in 
equations (45) and (46) will minimize the mean-squared error in c, without 
any regard to c, . Therefore, from equations (45)-(48), we can easiiy see that 
the MMSE interference canceller for signals from terminal unit 1 0 will 
consist of two different sets of weights a , and a, for c 2 and c„ respectively 
The weights for decoding signals from terminal 30 can be obtained in a 
similar fashion, as expected. Thus, the decoding of signals from terminal 
15 units 10 and 30 can be performed with a single subroutine MMSE.DECODE 
in decoder 25 as follows: 

(c,A r ) = MMSE.DECODE(r 1 ,r,,H 1 ,H ; ,G„G,,r) 
{ 

? = K r?] T 
M = HH* +— I 

■ r 

h i=[ h u hj,] = first column of H 
h, =[h, r 2 hj,] r = second column of H 
a' = M H h, , a; =M-'h 2 
c/ = argmin ja^-cj 2 , c, = argmin Jalr -c, 

A r =I a i ? -^ir+|a'r-c 2 " 2 
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With such a subroutine, both s and c can be estimated, as follows: 

( c,A)=MMSE.DECODE(r, ,r 2 ,H, ,H 2 ,G, ,G 2 ,D (49) 
(s,A)=MMSE.DECODE(r2,ri,G,,G 2 ,H u H2,r) (50) 

Similar to the zero-forcing case, simulation results reveal that the 
performance of the disclosed technique MMSE.DECODE is equivalent to 
that when only one terminal unit exists and the base station uses a single 
receive antenna which is equivalent to the performance of two branch MRRC 
diversity. However, this technique is also able to support two co-channel 
terminal units. In addition, when the SIR (signal-to-interference ratio, which 
is a ratio between the desired terminal power to the interfering terminal 
power) increases, the MMSE approach will have a better performance as 
15 compared to the ZF case (the ZF performance is almost the same as the 
performance of the MMSE approach at 0 dB SIR). 

Two-Step Interference Cancellation: Actually, additional 
improvement can be realized by employing a two-step interference 
cancellation approach using either the zero-forcing or the MMSE 
20 interference cancellation techniques disclosed above. Below, we will 
describe this approach based on the MMSE technique. However, as one 
might expect there a similar approach based on the zero-forcing technique. 
In this two-step approach, the receiver decodes signals from both terminals 
using the subroutine MMSE.DECODE disclosed above. Assuming that 
25 symbols from the terminal unit 10, c,, have been decoded correctly, the 
receiver can, then, perfectly cancel the contribution of the terminal unit 10 in 
the received signal vectors r, and r 2 . The receiver then uses x,. and x 2 , the 
received signal vectors after canceling signals from terminal unit 10, to re- 
decode symbols from terminal unit 30 s 0 using the optimum ML decoding 
30 rule in equation (10). Assuming that the symbols from terminal unit 10 has 
been decoded correctly, we can easily see that, the performance for terminal 
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unit 30 will be equivalent to that with 2 transmit and 2 receive antennas 
(which is equivalent to 4 branch MRC diversity). If we let 
A D = A r + A s denotes the overall uncertainty for c 6 and s 0 . The receiver 
then repeats the above step assuming that symbols from terminal unit 30 s, 
5 have been decoded correctly using the MMSE.DECODE subroutine. As 
before, the receiver cancels the contribution of terminal unit 30 in the 
received signal vectors r, and uses y, and y 2 , the received signal vectors 
after cancelling signals from terminal unit 30, to re-decode symbols from 
terminal unit 10 c, using the optimum ML decoding rule in equation (10). 

10 As before, assuming that symbols from terminal unit 30, the performance for 
terminal unit 10 will be equivalent to that with 2 transmit and 2 receive 
antennas. Similarly, let A ( = A r( + A J( denotes the overall uncertainty for c, 
and s r The receiver then compares the overall uncertainty and chooses the 
pair (c 0 ,s 0 ) if A 0 <A, and (c,,Sj) otherwise. The two-step interference 

15 cancellation approach based on the MMSE technqiue disclosed above is 
presented below in pseudo-code subroutine ILMMSE.DECODE. As we 
mentioned earlier, the techniques can be also used with the zero forcing 
approach. Below, we also present the pseudo code subroutine 
ILZF.DECODE for the two-step interference cancellation based on the zero- 

20 forcing approach. 
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(£,§)= II. MMSE.DECODE(r 1 ,r,,H,,H,,G 1 ,G 2 ,r) 

{ 

(c 0 ,A ro ) = MMSE.DECODE^.rj.H, ,H,,G,,G,,T) 
.X.-r.-H.rC.. , x 2 =r 2 -H 2 c„ 
s 0 = argmin (||x, -G, -s|| 2 +||x, -G 2 -s|| 2 ) 

seS 

A^^I^-G.-sf+K-G.-sf 
(Si , A Sil ) = MMSE. DECODE(r 2 ,r,,G 2 ,G,,H 1 ,H 2 ,r) 
y^r.-G.-s, , y 2 =r 2 -G 2 -s, 

c, = argmin (||y, -H, -cf +||y 2 -H 2 -cf ) 
A rJ =||y l -H 1 -cf + |y 2 -H 2 .cf 

If (A f . 0 + A J . 0 )<(A ( . 1 + A J . 1 ) => (c,s) = (£„,§,) 
Else (c,s) = (c,,s,) 

} 

(c,s)= II. ZF.DECODE(r l ,r 2 ,H„H 2 ,G„G 2 ) 

(c 0 ,A ro ) = ZF.DECODE(r„r 2 ,H 1 ,H 2 ,G 1 ,G 2 ) 
x,=r,-H, c 0 , x 2 = r,-H 2 c D 
s a = argmin (flx.-G, -sf +||x 2 -G 2 -sf ) . 

A^HK-G.-sf + IK-G.-sf 
(§„ A,.,) = ZF.DECODE(r,,r,,G 2 ,G,,H 2 ,H 1 ) 
y. = r i-G,-s 1 , y, = r, -G 2 -§, 

c, = argmin (||y, - H, • cf + ||y, - H, • c|f ) 
A^.^y.-H.-cf +||y J -H 2 .c|I I 

If (A,. <) +A J . 0 )<(A r , + A,,) (c,s) = (c 0 ,sj 
Else (c,s) = (£,,£,) 
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Interference Cancellation and ML Decoding: GENERAL CASE 
In the above basic case, we focused on the basic case where we 
assumed two co-channel terminals (K=2) each uses two transmit antennas 
(N=2). Both terminals communicate with a base station that is equipped with 
5 two transmit antennas (M=2). In this section we will consider the more 
general case of K > 2 co-channel terminals each is equipped with N > 2 
transmitting antennas, both terminals communicate with a base that has 
receive M > K antennas. We will develop similar interference cancellation 
and ML decoding scheme for this case. 
10 In a paper submitted to IEEE Transactions on Information Theory, 

Vahid Tarokh et al. extended the above space-time block coding scheme to 
the case when more than two antennas are used for transmission (N > 2). 
There, a technique for constructing space-time block codes (with similar 
properties to the simple scheme described above) was developed. It was also 
15 shown that for real constellations space-time block codes with transmission 
rate 1 can be constructed. However, for a general complex constellation the 
rate of transmission for these codes will be less than 1. 

In general, let us assume that the input information symbols to the 
transmitter are grouped into groups of Q symbols c l ,c 2i -^c Q . A space-time 
20 block code in this case will map the symbols c ii c 2y --,c Q into an N xL array 
C whose entries are made ± c x ,± c 2 , • * • ,± c Q and ± c\ ,± c* , • - ,± c Q . At time /, 
where 1< t < L, the r-th column of C is transmitted from the N antennas. In 
this case, the transmission rate of such code will be QIL. In general, for a 
rate QIL space-time block code (as constructed by V. Tarokh et al.) designed 
25 for transmit antenna, let r i9 r 2y - ,r L be the received signals at time 
t — 1,2,- ■ - 9 L . As before, we define the received signal vector as 

where the Lx 1 vector r can be written as 

r = H-c + t] 

30 and H is the LxQ channel matrix whose entries are 
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from±h^±h 2 ,'--,±h Ny ±h*,±hl,-",±h* N , and it is an orthogonal matrix, 
c = [c, c 2 •" c 0 ] r , and r\ is an Lxl zero-mean complex Gaussian random 

vector with covariance N Q • I which models the noise. The ML decoder in 
this case is similar to that in equation (6), that is 

5 c = arg mm||r - H • c|| 2 (53) 

and the uncertainty, A f , of the decoded symbols c is given by 

A r =||r-H-c| J . (54) 

As before, since the channel matrix H is orthogonal, the decoding rule in 
(53) decomposes into Q separate decoding rules for c n c 29 ~' 9 c Q . For 
10 example, assuming that the terminal unit uses 4 transmit antenna, a rate 4/8 
(i.e. it is a rate 14) space-time block code is given by 



V 








-c 3 


~ C A 


* 


m 


• 


C 2 


— > 


C 2 


C X 


C A 




• 
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Ca 
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C 3 




c y 


~ C A 


C l. 


c, . 


c 3 






- C 4- 






C i 


-c. 











In this case, at time r=l c l9 c 2r c Z9 c A are transmitted from antenna 1 through 4, 
respectively. At time r=2, -c 2 , c } -c A , c 3 are transmitted from antenna 1 
15 through 4, respectively, and so on. For this example, let ^,r 2 ,--*,r g be the 
received signals at time r = l,2,---,8. Define the received signal vector 

r = [r x r 2 r 3 r 4 r* r* r* r 8 *] r . In this case, we can write the received 
signal vector r can be written as 

r = H c + r| (55) 
20 where r| is the 8x1 AWGN noise vector and H is the 8 x 4 channel matrix 
given by: 
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-K 
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hi 
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-K 
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-K 


-hi 


hi 


K 


hi 


-hi 


-k 



(56) 



We can immediately notice that the matrix H is orthogonal, that is 
H'H = D h • I where D h = ^|^,|~ and I is a 4x4 identity matrix. 

' /-"I ' 

5 

Let us now assume a multi-user environment with K co-channel 
synchronized terminals. Each terminal uses a rate QIL space-time block code 
with N transmit antenna (as constructed by V. Tarokh et al). The base station 
uses M> K antennas for reception. We can write the received signal vector 
1 0 at the m-th receive antenna as 

K 

r m =XH, m -c, + ti„ 1 , m = \X---,M (57) 

where is the LxQ Avth user channel matrix to antenna m , 

c k ~\c k} c k2 c kQ ] T is the Qxl information symbols vector for k-th 

user, and T\ m is the Lxl noise vector. The entries of the k-th user channel 
15 matrix are from ±h kM ^±h k ^-,±h k ^ N and ±hl m]9 ±hl m2 r''+hl m „ J 
where h k m n is the complex channel gain between transmit antenna n of the k- 
th user and receive antenna m. As we stated earlier, the matrix is 
orthogonal. 

Zero-Forcing IC and ML Decoding : Without loss of generality, let us 
20 assume that we are interested in suppressing signals from co-channel 
terminals 2,3,---,AT while decoding signals from the first terminal. This can 
be done in a successive manner as follows. 
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First, let us define r/ H 0) = r m . Let us assume that we start by canceling 
out the contributions of the K~th terminal. We can use the M-th antenna 
received signal vector r u to cancel out the contribution of the K-th terminal 
in the remaining AM received signal vectors by forming the modified 
5 received signal vectors r^\m = 1, --,A/- 1 as follows: 

where is the generalized inverse of the channel matrix Yi hn and is given 
by 

10 We can easily verify that Hj ffl H Am = I, where I is the QxQ identity matrix. 
We can easily verify that the modified received signal vectors 

m = l,-*-,M — 1, do not contain any signal contribution due to the AT-th < ; 
user. Moreover, we can easily verify that can be written as 

li'^SH^-c^ + Ti^.. m = l,2,..-,M-l (60) 

15 where H^and -q^ are given by 

H^ = H^-H^(H^) + H^, m = l,2,-,M-l (61) 

^ = ^?-^ 0 l{nZ) + <\ ™ = 1,2,-,M-1 (62) 
Moreover, it can be shown that for codes constructed by V. Tarokh et al, the 
modified channel matrix h>/// have exactly the same structure as that of 
20 . That is, the entries of the k-th user modified channel matrix are 

from ±C,l±*ii 2 r-s±AjI^ and ±^^^. 2 r--.±A^. where is 
the modified complex channel gain between transmit antenna n of the /r-th 
user and receive antenna m, m=l,„.,M-l. Hence, the modified channel 
matrix H ( A 2 will be orthogonal as well. 
25 It can then be observed that the expression for the M-l modified 

received signal vector r£ } in equation (60) is the same as that in equation 
(57) except that we now have one less interfering terminal. In a similar 
fashion, we can cancel out the contributions of terminal K-\ and obtain M-2 
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modified received signal vectors r { J\m = l, --,M-2 that do not contain any 
contributions from terminals AT-th and AM. In general, after stage j\ where 
y=l,2,..,AM contributions from terminals K^K-ly.^JC-j+l are canceled out 
and we are left with M-j modified received signal vectors 
5 r^,m = \,-*,M — j\j : = 1,2,- -,AT-1 that contain signals due to terminals 
1 ,2, . . . ,K-j only. In this case, we will have 

ri/^XH^-c.+Ti^, m = l,2,...,M-y (63) 

where H ( A ;^ and rf^ are given by 

r w = r u-o _ H o-n (H (y : ., A// )♦ r o_- ; , m = 1>2i . . . M _ j (M) 

= H ( £'» - H^H£^ + H<^ , 1 < m < M-j, \<k<K-j (65) 

^ = rir' -^LW^-jf^j, rn = l,2,-,M-7 (66) 
This process is repeated until we are left with M-K+l modified 
15 received signal vectors r { *~ x \m = A*- AT + 1 that contain only 
contributions due to the first terminal. In this case we will have 

r^-'^HiJ^-^+.Ti^ m»l,2,..., M-K + l (67) 
which contains signals due to the first terminal only. Similarly, the modified 
channel matrix H** -0 , rn = 1,2,- • M - K + 1 , will have a similar structure and 
20 is also orthogonal. Hence, it is straight forward to see that the ML decoder 
for signals from the first terminal is given by 

c, = argmin"x ' ft™ - H£-'» • c| (68) 

and the corresponding uncertainty will be given by . . . 

X llnT'-H^'Hl (69) 

25 Similarly, since the modified channel matrices H|*~'\ 1 < m < M-K + l are 
orthogonal, as before, the decoding rule in (68) will decompose into Q 
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separate rules for decoding c,,,c 12 ,* ~<,c XQ . We may observe that the basic 
case for zero-forcing IC and ML decoding that we discussed in detail 
earlier is a special case of the above approach. 

(c, A) = G_ ZFDECODE({r„, } . {H lm }^ K I ^ M ) 
{ 

r l m 0) = r m , \<m<M 

H< 0 i=H,, m , \<m<M, \<k<K 

for y = l-> K-\ 

M,= M-j,K t = K-j 

for / = 1 — > M j 

H\f] = H^'> - H^(H^ / ) + H^> , \<k< Kj 

end 

end 

c = argmin Af £ , ||rr'>-H < | -'>£|r 

»i=l 

} 

The above-disclosed technique can be easily implemented within a 
detector 25 that comprises a stored program general purpose processor. 
Specifically, a subroutine (c,A)=G,ZFDECODE({r,„} ^ tSAf ,{H*, M } 
\<k<K,\< tn <M ) can be installed which returns the values c,A in response to 
10 submitted inputs {r,„} u,„<^ and {H km } \<k<x,\<m<M , as shown above. 

Minimum Mean-Squared Error IC and ML Decoding Scheme: The 
MMSE IC and ML decoding in the general case can be developed in a 
similar fashion as follows. We recall the received signal vector at the m-th 
15 receive antenna in equation (57) 
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10 



15 



•",„ = Xh,,„-c, +ti„„ m = \,2,--,M 

This can be written in a matrix form as in equation (38) 

r== H c + n 

where r = [r t r r 2 r r^j is a MZxl vector, c=[cf c[ 

QKx 1 a vector, n = [*q J r\l * • * t^,] is a MZ,x 1 vector, and 



z T K ] T \s 



H I2 H 22 



H 



K2 



H, 



is the MLx OK channel matrix. As before, we redefine the vector r as 
r — [ r i r 2 "' r Mif • As before, we assume that we are interested in 
decoding the symbols transmitted from terminal 1 c x x , c, 2 , • - • , c, Q . As before, 
when seeking to detect and decode signals c,,,c I2 ,*--,c ig by minimizing a 
mean-squared error criterion, the goal is find a linear combination of the 
received signals such that the mean-squared error in detecting the signals 
c U9 c l29 -" 9 c } Q is minimized. In general terms, this can be expressed by an 
error cost function that is to be minimized, such as the function 

Aa.fi) = 

<=l 7=1 

Similarly, as before we can see that one of the fi p \ < j <Q must be set to 1 
or else we get an all zero solution for a and fi. Consider the case where we 
set fij = 1 . Hence, in this case, the criteria to be minimized is 



~|a.r-l-e,r 



-foXjTj-cJ , \ <j<Q 



(70) 



(71) 



(72) 



(73) 



(74) 



20 where 



a j = [«,p a j , t2 , ■ • ■ , a jLM , • • • -/?._, -J3 /+I ,~ -J3 Q ] = [a . -fi, ] 



(75) 
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7 J=[ r J c i 



1 1 




(76) 



" 10 



15 
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If we follow the same steps as in the basic case, we arrive at the conclusion 



where h ; is the >th column of the channel matrix H, and M = HH* + — I, is 

an MLxML matrix, T is the signal to noise ratio, and I is the MLxML 
identity matrix. 

In this case, as before, the error in decoding the y-th symbol Cjy will be 
minimized without any regard to the other symbols. Hence, the MMSE-IC 
and ML decoder will consist of Q different combiners, one for each symbol. 
It should be clear now that the MMSI-IC solution for the general case is a 
straight forward extension to the basic case shown earlier. The MMSE-IC 
solution for the general case can be implemented using the subroutine 
G MMSE.DECODE shown below. 



that 



= 1 i = y 
a] =M-'h., \ <j<Q 



(77) 



(78) 
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(c,A) = G.MMSEDECODE({r m }^ A/ ,{H^} 



•0 



H = 



H ll H 2I 

^12 ^22 ^K2 
.^IW ^2 A/ * H/CW. 

M « HH*+ — I 

r 

for 7 = l-*Q 

h . : == 7-th column of H 

a*=M~'h.. 



Cj = argmm a 7 r-.c. 



end 

c = [c, c 2 

Q 

7-1 
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We claim: 

1. A method for decoding a set of M signals received at an input 
interface from a plurality of terminal units that transmit on a given channel, 
comprising the steps of: 

processing said M signals, received over L time intervals, where L is 
an integer, with signals having components related to channel coefficients 
between transmitting points of said terminal units and an input interface, to 
detect signals transmitted from each of said terminal units by canceling 
interference from K of said terminal units, where M and K are integers such 
that M> 2 and K < M , and to identity probable signals transmitted by said 
terminal units through maximum likelihood detection; and 

applying signals developed through said maximum likelihood 
detection to a location, from which the signals may be applied to post 
processing that culminates in signals adapted for delivery to users. 

2- The method of claim 1 where transmissions of said terminal units 
are synchronized. 

3. The method of claim 1 where each of the K terminal units employs 
at least two transmitting antennas, and said input interface includes M 
receiving antennas. 

4. The method of claim 1 where each of the K terminal units employs 
N transmitting antennas, where TV is greater than 1 , and said input interface 
includes M receiving antennas. 

5. The method of claim 4, where said interference canceling involves 
algebraic operations carried out on said M received signals with said signals 
having components related to channel coefficients between transmitting 
points of said terminal units, to form a plurality of signals, each of which is 



WO 99/18682 



PCT/US98/20982 



-28- 



substantially unrelated to signals sent by all but one of the terminal units. 



6. The method of claim 4 where said K—M=2 and 1 < N < 2 , 



7. The method of claim 6 where said processing comprises 
multiplying said vector of received signals by conditioning matrix 



-g.g;' 



i 

-h 2 h;' i 

*21 ^22 
Jh2 ~ A 2.. 



, where H. = 



! .. K. 1 

i2 — VlJ 



G,= 



H, = 



$\2 S\ I ^ 

, I is the diagonal matrix, h i} is a 



§2 1 8 22 
* ♦ 

channel coefficient estimate between a transmitting antenna / of a first 
10 transmitting unit and a recei ve antenna j of said M receiving antennas, and 
g u is a channel coefficient estimate between a transmitting antenna / of a 

second transmitting unit and a receive antennay of said M receiving 
antennas. 



15 



8. The method of claim 7 where L=2. 



20 



9. The method of claim 7 where said maximum likelihood detection 
minimizes a metric |jr; — S - c[| 2 over all potential code vectors c 9 where r; is 
one of said plurality of signals formed by said step of multiplying, and 
H = H, — GjG-j'H, . 



10. The method of claim 9 where said maximum likelihood detection 
develops an uncertainty measure. 



25 



11. The method of claim 10 where the uncertainty measure is given 



byA r = ||l;-H.c" 2 



12. The method of claim 6 where said processing executes a 
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subroutine ZF. DECODE(r t , r, , U x , U y , G , , G , ) which returns an output c 
by executing the following calculations 

H = H r -G r G;H v) and 
c = argmin llr -H-cIf, 

where 

is the first argument of the ZF.DECODE subroutine and corresponds to 
a vector of signals received over said L time intervals at receiving 
antenna jc, 

r v is the second argument of the ZF.DECODE subroutine and corresponds 
to a vector of signals received over said L time intervals at receiving 
antenna y y 

H x is the third argument of the ZF.DECODE subroutine and corresponds to 
a matrix comprising channel coefficients between transmitting antenna* 
of a first terminal unit and said M receiving antennas, 
is the fourth argument of the ZF.DECODE subroutine and corresponds 
to a matrix comprising channel coefficients between transmitting 
antenna of said first terminal unit and said M receiving antennas, 

G r is the fifth argument of the ZF.DECODE subroutine and corresponds to 
a matrix comprising channel coefficients between transmitting antenna x 
of a second terminal unit and said M receiving antennas, 

G y is the sixth argument of the ZF.DECODE subroutine and corresponds to 
a matrix comprising channel coefficients between transmitting antenna^ 
of said second terminal unit and said M receiving antennas, and 

c is a selected code vector c from a set of code vectors C as the code 

transmitted by the terminal related to the second and third arguments of the 

ZF.DECODE subroutine. 



13. The method of claim 12 where said processing comprises 
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executing the ZF.DECODE subroutine a first time with argument 
(r^Fj.H^H^G^Gj), and executing the ZF.DECODE subroutine a second 
time with argument (r 2) r 1 ,H 2 ,H I ,G 2 ,G 1 ). 

5 14. The method of claim 12 where said subroutine ZF.DECODE also 

returns an uncertainty measure output A c by carrying out the calculation 

A c =[r-H.f . 

15. The method of claim 14 where said processing comprises the 
10 steps of: 

executing the ZF.DECODE subroutine a first time, with argument 
(r l ,r 2t H l ,H 29 G |t G 2 ), to obtain a first version of a probable vector signal of 
a first terminal unit, c 0 , and a measure of uncertainty, A r 0 ; 

obtaining a first version of a probable vector signal of a second 
15 terminal unit, s 0 , and a measure of uncertainty A s 0 ; 

executing the ZF.DECODE subroutine a second time, with argument 
(r 2 ,r,,G 2 ,G p H 3 ,H,), to obtain a second version of a probable vector signal 
of the second terminal unit, s G , and a measure of uncertainty, A s , ; 

obtaining a second version of a probable vector signal of a first 
20 terminal unit, c 0 , and a measure of uncertainty A c , ; 

selecting said first version of said probable vector signals if 
(A r 0 + A J 0 ) < (A c , + A 5 ,) , and selecting said second version of said 
probable vector signal based otherwise. 

25 16. The method of claim 1 where M > 2 and K > 2, and where said 

processing comprises the steps of: 

combining said M received signals for form a signal that is 
essentially independent of transmissions of all but one of the terminal units, 
performing a maximum likelihood detection to identify a probably 
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signal vector transmitted by said one of the terminal units, 

modifying remaining ones of said M received signals to account for 

contribution of the probable signals vector, and 

returning to said step of combining as long as there is a terminal unit 
5 whose transmission is to be detected, to select another one of said terminal 

units whose signal is to be detected with said step of maximum likelihood 

detection. 

17. The method of claim 1 where M > 2 and K > 2, and where said 
10 processing comprises executing a subroutine 

(c,A) = G ZFDECODE({r }, 4 ,,{H, ), „) 
{ ' 
r™=r H , \<m^M 

= H ljn , \Sm%M, iSk&K 
for j = K-\ 

M J = M-j,K j =K-j 
for i = l—> Mj 

r/>'=^-''-HV;;'>(H^»j + r<;;'» 

end 

end 

M-K+) 

c = argmin £ |^*-'» -B\?> -cf 

C|€C m=l 

A=Y'lk: A: " ,) -H!r , -ci 2 

m-l 

} 

where 

rj? } is a received signal at antenna /«, as modified in step (a); 

is a matrix of channel coefficients between terminal unit k and receive 
15 antenna m, at step a; and 

(Hi a j,) + is the generalized inverse of . 
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18. The method of claim 1 where said processing form a vector of 
signals from said M signals and said processing includes pre-multiplying said 
vector by a preconditioning matrix having components related to said 
channel coefficients. 

19. The method of claim 18 where said components of the 
conditioning matrix are related to estimates of said channel coefficients. 

20. A detector responsive to M received signals, where M>2, from a 
plurality of terminal units and canceling interference from K of said terminal 
units, where K < M i comprising: 

a processor for transforming said M received signals through 
algebraic operations involving components related to channel coefficients 
between transmitting antennas of said terminal units and said M received 
signal to form a plurality of signals, to each of which is substantially 
unrelated to signals sent by all but one of the terminal units, where M, N, and 
K are positive integers, and to perform maximum likelihood detection 
calculations on said plurality of signals; and 

a memory coupled to said processor for storing information for 
controlling said processor. 

21. A detector responsive to M received signals, where M >2, from a 
plurality of terminal units and canceling interference from K of said terminal 
units, where K < M , comprising: 

a first means responsive to said M received signals, for transforming 
said M received signals through algebraic operations involving components 
related to channel coefficients between transmitting antennas of said terminal 
units and said M received signal to form a plurality of signals, each of which 
is substantially unrelated to signals sent by all but one of the terminal units, 
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where M, N, and K are positive integers; 

a second means, responsive to said first means, for performing 
maximum likelihood detection calculations on said plurality of signals; and 

means for applying information developed by said second means to a 
location accessible for further processing and delivery to a user. 

22. The detector of claim 20 where the processor also calculates an 
uncertainty measure associated with each maximum likelihood detection. 

23. The detector of claim 20 where said processor forms said 
plurality of signals one at a time, and with each formed signal of said 
plurality of signals the processor selects a probable transmitted signal vector 
of one of said terminals and evaluates an uncertainty measure. 

24. The detector of claim 20 where said processor, in the course of 
forming one of said signals from said M received signals, substantially 
nullifies contribution of transmissions from terminals for which a probably 
transmitted signal vector has been selected. 

25. The detector of claim 20 where said processor: 

forms a first version of a first of said plurality of signals, selects a 
first version of a probable transmitted signal vector of a first of said 
terminals, and evaluates a first uncertainty measure, then 

forms a first version of a second of said plurality of signals and in the 
course of said forming said first version of said second of said plurality of 
signals nullifies contribution of transmissions of said first version of a 
probable transmitted signal vector of said first of said terminals, selects a 
first version of a probable transmitted signal vector of a second of said 
terminals, and evaluates a second uncertainty measure, then 

forms a second version of said second of said plurality of signals, 
without said nullifying, selects a second version of a probable transmitted 
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signal vector of said second of said terminals, and evaluates a third 
uncertainty measure, then 

forms a second version of said first of said plurality of signals and in 
the course of said forming said second version of said first of said plurality 
of signals nullifies contribution of transmissions of said second version of a 
probable transmitted signal vector of said second terminal, selects a second 
version of a probable transmitted signal vector of said first of said terminals, 
and evaluates a fourth uncertainty measure, then 

forms a first combination involving said first and second uncertainty 
measures, forms a second combination involving said third and fourth 
uncertainty measures, and selects either said first version or said second 
version of the probable transmitted signal vector of said first of said 
terminals and the probable transmitted signal vector of said second of said 
terminals based on whether said first combination is greater than said second 
combination. 

26. The detector of claim 20 where M=N=K=2. 

27. The detector of claim 20 where M>2, and A>2. 

28. The detector of claim 20 where M>2 y K>2, and 7V>229. The 
detector of claim 20 where the processor also estimates parameters of 
channel through which said M signals traverse. 

30. The detector of claim 22 further comprising M antennas and 
amplification and demodulation circuitry interposed between said M 
antennas and said processor. 

31. The detector of claim 20 further comprising a subroutine stored 
in said memory and executed by said processor, which subroutine performs 
said transforming of said M received signals. 
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32. The detector of claim 24 where each execution of said subroutine 
develops one of said signals that is substantially unrelated to signals sent by 
5 all but one of the terminal units. 

33 The detector of claim 25 where said subroutine is executed not 
less than K times. 

10 34, The detector of claim 20 further comprising a subroutine stored 

in said memory and executed by said processor, which subroutine performs 
said transforming of said M received signals, followed by a maximum 
likelihood detection. 

15 35. The detector of claim 20 further comprising means for estimating 

channel parameters. 

36. The detector of claim 20 where said processor forms said 
plurality of signals by forming a vector from said M received signals and 
20 pre-multiplying said vector by a conditioning matrix of coefficients that are 
related to channel coefficients between said transmitting antennas and an 
input interface.. 



25 



37. The detector of claim 13 where said conditioning matrix is 



[-HW 



i 



, where H, = 



A*2 -*fu 



8 

.8*2 



1 1 Sl2 I 

12 Sll J 





"#2. 


§22 ' 








.£2 





, I is the diagonal matrix, h u is a 



channel coefficient estimate between a transmitting antenna i of a first 
transmitting unit and a receive antenna j of said M receiving antennas, and 
g r is a channel coefficient estimate between a transmitting antenna / of a 
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second transmitting unit and a receive antenna j of said M receiving 
antennas. 

38. The detector of claim 37 where said processor also performs 
5 maximum likelihood detection on said plurality of signals. 

39. The detector of claim 38, where said input interface comprises M 
receiving antennas. 
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